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Recent work by Teifel and Mahler [Eur. Phys. J. B 75, 275 (2010)] raises le- 
gitimate concerns regarding the validity of quantum nonequilibrium work relations 
in processes involving moving hard walls. We study this issue in the context of 
the rapidly expanding one-dimensional quantum piston. Utilizing exact solutions 
of the time-dependent Schrodinger equation, we find that the evolution of the wave 
function can be decomposed into static and dynamic components, which have sim- 
ple semiclassical interpretations in terms of particle-piston collisions. We show that 
nonequilibrium work relations remains valid at any finite piston speed, provided 
both components are included, and we study explicitly the work distribution for this 
model system. 

PACS numbers: 05.70.Ln, 05.30.-d, 05.40.-a. 05.90.+m. 

In the past two decades, much attention has been devoted to theoretical predictions and 
experimental investigations regarding the fluctuations of small systems away from thermal 
equilibrium. These predictions include the nonequilbrium work relation [H |2] 




(1) 



and the corresponding fluctuation theorem derived by Crooks [3H5] 




(2) 
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which pertain to the work (1^) performed on a system driven out of equihbrium. (See Ref. [6] 
for details and a recent review of these and related results.) Most of the research in this 
area has concerned systems evolving under classical deterministic or stochastic dynamics. 
However, the past few years have seen increased focus on the fluctuations of quantum systems 
driven away from equilibrium [7]. 

While the derivation of Eq. [T]for an isolated quantum system is straightforward and rests 
on familiar properties of unitary evolution [HHTT] , Teifel and Mahler (TM) [12] have recently 
presented a calculation suggesting that Eq. [l] (and by extension, Eq. [2]) might be violated 
for the one-dimensional quantum piston. In this familiar model system, the wavefunction 
describing a particle inside a box evolves in time as the length of the box is increased (Fig. 
[T| or decreased. Although TM focus specifically on this simple model, their analysis has 
broader implications, raising the possibility that Eqs. [T] and |2] might generically be violated 
for processes involving the motion of hard walls. In such situations the system's Hilbert 
space changes with time, and questions of unitarity must be handled with care. This feature 
has a classical counterpart, emphasized by Sung [13]: the phase space accessible to a classical 
particle confined by hard walls changes with time as those walls move. 

In the classical setting, processes involving moving boundaries have proven to be instruc- 
tive [T3HT7] . deepening our understanding of nonequilibrium work relations by highlighting 
apparent paradoxes and counterintuitive features. In this paper we use exact solutions of 
the time-dependent Schrodinger equation [TH] to investigate the validity of Eq. [l] for an 
expanding quantum piston. 

In what follows, we first sketch the usual derivation of Eq. [T] for an isolated quantum 
system (Eqs. |3]-[8]), as well as an apparent counter-argument which suggests that Eq. [l]is 
violated for the quantum piston (Eqs. [9|-[lT|. We then apply the exact results of Ref. [18] to 
the case in which the piston moves outward at speed v. We find that Eq. [T]is valid for any 
finite pulling speed, which seems to contradict the analysis in Eqs.|9]- 11 We then consider 
the limit f — )■ oo, and we find that the apparent discrepancy has an appealing semiclassical 
interpretation that parallels the purely classical analyses of Refs. [T314T6] . 

Consider a quantum system whose parameter-dependent Hamiltonian H'^ has a discrete 
energy spectrum: 

H^\m^) = E^\m^), (3) 
with m = 0, 1, 2 ■ ■ ■ . We use superscripts to indicate the value of the externally controlled 
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A = A 




A = B 



FIG. 1: Schematic depiction of a quantum piston. A quantum particle is confined by hard walls, 
one of which acts as an externally controlled piston. We focus on the case in which the piston is 
pulled outward at a speed that is much greater than the initial thermal speed of the particle. 

parameter, A, which for the case of the quantum piston is the position of the piston itself, 
equivalently the length of the box. Now imagine that this system is subjected to the following 
process. (1) With the parameter fixed at X = A, the system is equilibrated with a reservoir 
at temperature the reservoir is then disconnected and the energy of the system is 
measured. At this point the system is in a pure state \m^), set by the outcome of the energy 
measurement. (2) The system now evolves under Schrodinger's equation, from time t = to 
t = r, as the parameter is varied from A to B according to a schedule, or protocol, Xt- The 
energy is then measured once more, resulting in "collapse" into an eigenstate \n^) of . 
Following Refs. [814TT] we identify the work performed on the system with the net change in 
its energy: 

W^E^- Et (4) 

By repeating this process, we generate an ensemble of realizations, each defined by an initial 
state \m^) and a final state \n^). The initial states are distributed according to 

Pr^ = ^e-^^^, (5) 

where Za = Yl,m. e"^^™ = e~^^^ is the partition function, and the final states according to 
the conditional distribution 
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where tj is the time-evolution operator that describes evolution under Schrodinger's equation 
from t = to r. Combining Eqs. |4]-|6| the left side of Eq. [T]can now be evaluated: 



m n n 

where 



5r), — 



At this point, one normally argues that the sum \m^) (m^| is the identity operator, 
hence s„ = 1 and the right side of Eq. [t] becomes Zb/Za = e~^^^ , completing the proof. 

Teifel and Mahler [12] correctly point out that this argument requires care if the eigen- 
states of do not span the Hilbert space of . For a quantum piston whose length is 
increased from Aq = A to = -B at speed f , the states \m^) are restricted to the interval 
< X < A, whereas the final Hilbert space supports states extending over the wider interval 
< a; < -B. If il){x) = (xlip) is a wave function belonging to the Hilbert space of H^, then 
the operator Ylm \^^) {^^\ effectively "chops off" a portion of this wavefunction: 

^(a;|m^)(m^|^) = e(v4-x)^(x), (9) 

m 

where 9{-) is the unit step function. We conclude that Ylm 1^^) {^^\ is not the identity 
operator when it acts in the Hilbert space spanned by eigenstates of . Hence the derivation 
described in the previous paragraph does not automatically apply to the quantum piston, 
and this raises concerns regarding the validity of Eq. [T]in that context. 

As a limiting case, let us analyze the infinitely fast expansion of the piston, f — )■ oo. The 
sudden approximation [19] suggests that the wave function then remains in its initial state, 

lim U\m^) = |m^). (10) 

v—^oo 

Combining Eqs. |8 10 leads to 



lims„ = y(n^|m^)(m^|n^)= Tdx = - - ^^^^^^^^ < 1, (11) 

m=l •'^^ 



where r = BjA and the wavefunction 



0„(x;A) = ^-sm(— j (12) 
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describes the n'th eigenstate of H^. (The notation = indicates that we question the vahdity 



of the first step in Eq. 11 ) Substitution of Eq. 11 (s„ < 1) into Eq. [Tjimphes a violation 
of Eq. [1} In the opposite hmit, namely adiahatic expansion, f — )■ 0, TM find that Eq. [T]is 
satisfied. These considerations suggest that for the expansion of a quantum piston at finite 
speed f , Eq. [T] is only approximately valid, but the approximation becomes exact in the 
adiabatic limit, f — )■ 0. 

In what follows we will argue that in fact s„ = 1 for all finite values of n and v, and 
therefore 

lim Sn = 1, (13) 

v—^oo 

in contradiction with Eq. VL By Eq. [7], our conclusion implies that Eq. [T] is valid for any 
finite piston speed. 

For a quantum piston expanding at speed v from an initial length Aq = ^, a set of 
independent solutions to the time-dependent Schrodinger equation can be written as ^18] 

1 



^i{x,t) = exp 



-Mvx' - Ef'At 



<Pi{x] At), 



/ = 1,2, 



(14) 



where M denotes the mass of the particle, and Ef = fn'^h'^/2MA'^ is the Uth eigenenergy 
of the system at t = 0. The wavefunctions ^i{x,t) form a complete orthonormal set, 
= 6ki, but are not eigenstates of H^K (The (pi's defined in Eq. 12 are the eigenstates.) 
A general solution of the time-dependent Schrodinger equation takes the form 

oo 

^{x,t) = Y,ciMx,t), (15) 



1=1 

where the time-independent coefficients q are set by the initial wave function: 

"A 



Cl 



<l>;(a;,0)^(a;,0) dx. 



(16) 



For initial conditions |^'(0)) = |m^) these coefficients are (setting h = M = 1) 



ci{m) 



A Jo 



lirx 
~A 



sm 



i"rmTX\ 



\ A 



-) 



dx, 



(17a) 



and the transition matrix element to the state \n^) at the final time r is 



oo „B 

'|f/|m^) = (n-^l^(r)) = Vc^m) / ^^(a;; 5) $^(0;, r) dx. (17b) 

1=1 -^0 
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Eqs.[6]and 17 give the transition probability P{n^\m'^), in terms of one-dimensional integrals 
that are easily computed numerically. This transition probability satisfies normalization: 

2 







/ dx 


{x\U\m^) 


Jo 





1. 



(18) 



Although we have considered the expansion of a quantum piston, Eq. [14] is equally valid 
for compression [18] . By reversing the roles of A and B and the roles of m and n and by 



replacing v with —v (in Eq. 17) we obtain the transition probability P(m \n ) from the 
n'th eigenstate of to the m'th eigenstate of H^, where the notation P indicates the 
compression process. This transition probability also satisfies normalization: 

m 

where II' is the time-evolution operator for the compression process. 



.A 




\n^) = dx 


(x|[/'|n^) 


Jo 





(19) 



In the Appendix, we provide explicit expressions for P(n \m ) and P{m \n ), and using 



these expressions we directly verify the relation 



P{n^\m'') = P(m^|n^) 



(20) 



It should be clear that this relation is precisely what we expect from time-reversal invariance 
(f/' = f/t), see e.g. Eq. 56 of Ref. [7]. Using Eq. 
initial states in Eq. [8] into a sum over final states: 



20 



we can now transform the sum over 



Y^P{n^\m^) = ^P(m^ln^) = 1 



(21) 



using Eq. [19] in the last step. Since this result is independent of f , we conclude that Eq. [T] 
is valid at any finite speed of expansion. 

To obtain Eq. [2] by similar means, we follow Tasaki [0] and write explicit expressions for 
the forward and reverse work distributions (corresponding to piston expansion and compres- 
sion, respectively): 



m n 



m \n 



(22) 



For every realization m"^ — )■ that gives a particular work value during the forward process, 
there is a corresponding realization — )■ that gives the opposite work value during the 



reverse process. Combining this observation with Eqs. 20 and 22 we obtain Eq. ^ [9]. 
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Up to this point we have used the symmetry relation, Eq. 20, to show that s„ = 1 for 
any finite speed f , and therefore that Eqs. [T] and |2] remain vahd for the quantum piston. 



However, this analysis does not yet explain why Eq. 11 gives a contradictory result in the 
limit f — > oo. To address this issue, in the following paragraphs we present numerical 
evidence that the value of s„ is naturally expressed as the sum of a static and a dynamic 
contribution, refiected in the two-peak structure seen in Figs. [2(iO] -[2(^ The sum of these 



contributions is unity for any finite v (as per Eq. 21), but Eq. 11 accounts only for the 
static contribution, thus giving Sn'^ '^^ < 1. Here and in the following discussion, we use the 
notation Sn^ '^^ to denote the value for s„ predicted (incorrectly!) by Eq. [ll| in the limit 
f — 7- oo. After presenting the numerical results, we suggest a semiclassical interpretation in 
terms of piston-particle collisions. 



We have used Eq. 17 to evaluate P(n \m ) numerically. In Fig. ^ this quantity is plotted 



for fixed final state function of initial state m = 1, 2, ■ ■ ■ , for piston expansion 

from A = 1.0 to B = 2.0 at various speeds: v = 10, 100 and 500. The plot reveals a 
two-peak structure. The left peak, near m = 2, remains approximately independent of v, 
whereas the right peak is located near m = 2vA/tt; thus with increasing v the right peak 
shifts further rightward. (Note the change of scale in the plots.) We will refer to the left 
and right peaks as the static and dynamic components, respectively. We can decompose the 
value of Sn, with n = 3 in our case, into contributions from these components: 



J2 P{n^\m^), 



-.R 



J2 P{n^\m^). 



(23) 



m<m* 



m>m* 



Here m* is the value of m at which P(n^|m^) is minimized in the region between the two 
peaks. Table |I] lists the values of these contributions, obtained by numerical evaluation of 



the integrals in Eq. 17, as well as their sum, Sn- Note that Sn = 1.000 at all three speeds, in 



agreement with Eq. 21 



Let us now rewrite Eq. 11 



as 



A 2 
dx (pKx; B)(l)mix] A) 



(24) 



lim Sn = y^in^lm^) {m^\n^) = 

m m 

We can interpret the overlap 0(n-^|m"^) = | (n'^|m'^) as the probability to end in state 
\n^) after the measurement of the final energy, when starting from state |m^), under the 
assumption that the wave function remains unchanged during the sudden expansion. This 
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FIG. 2: P{n^\m^) is plotted as a function of m at fixed n = 3 (open circles), revealing a two-peak 
structure, with the left peak around m = 2 and the right peak near m = 2vA/7r. We refer to these 
peaks as the static and dynamic components, respectively. Also plotted is the quantity 0{n^\m^) 
(red points), which displays only a single peak around m = 2. Note that for v = 100 and v = 500 
the single peak of 0{n^\m^) is virtually identical to the static component of P(n^\m^). 



assumption amounts to a literal intepretation of the sudden approximation, Eq. 10 Using 



Eq. 12 to evaluate the integral, in Fig. [2] we have also plotted Oln^lm^), which exhibits 
a single peak around m = 2. We observe that the larger the value of v, the more closely 
0{n^\m'^) resembles the left peak of P{n^\m'^); indeed at f = 100 and 500 they are virtually 



identical. These empirical observations suggest that Eq. 11 captures only the contribution 



to Sn from the static component s^, while missing the contribution from the dynamic com- 
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V = IQ 


V = 100 


V = 500 


t> — )• oo 


Eq.fTT} 








0.500 




0.644 


0.499 


0.500 






0.356 


0.501 


0.500 




•Sn = •Sn + •Sn 


1.000 


1.000 


1.000 





TABLE I: Static (L) and dynamic (i?) contributions to Sn=3, as well as the asymptotic value of Sr, 
predicted by Eq. [TTj for piston expansion from ^ = 1.0 to = 2.0. 



ponent s^. 

Quantitatively, Sn*^ '^ = 0.5 for A = 1.0, B = 2.0, and n = 3, whereas the data in Tablejl] 
suggest that the static contribution approaches 0.5 as f — ?■ oo. Moreover, Table [n| hsts 
these quantities for the case A = 1.0, B = 1.485, and n = 3, with and again calculated 
using Eq. 17 Once again we find that s^ + s^ = 1.000 at all speeds, and — )■ Sn'^ 'I^^' ~ 0.667 



as f — )■ oo. These findings are consistent with our hypothesis that Eq. 11 refiects only the 
static and not the dynamic contribution to Sn- 





v = 10 


V = 100 


V = 500 


u — )■ oo 


Eq.([TT} 








0.667 


•'n 


0.638 


0.667 


0.667 




qR 


0.362 


0.333 


0.333 




•Sn •Sn ~l" •Sn 


1.000 


1.000 


1.000 





TABLE II: Same as Table |IJ but for expansion from A = 1.0 to i? = 1.485. 



We now build a semiclassical interpretation to reinforce these conclusions. For A = 1.0, 
B = 2.0, and piston speed v = 100, consider the value P(3^|64^), corresponding to the 



right peak in Fig. 2(b), This gives the probability to end in state |3'^), starting from state 
164"^), during the expansion process. Semiclassically, the initial state |$(0)) = |64^) can be 
imagined as a particle moving with speed 



\u\ 



2-E'm=64 



rnvr 



200 



(25) 



between two hard walls. At t = 0, when the piston begins to move rightward with speed 
V = 100, the particle is moving either leftward {u ~ —200) or rightward {u ~ +200), 
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with equal likelihood. In the latter case, the particle will collide once with the receding 
piston, losing approximately all of its kinetic energy. The final state |$(t)) will then contain 
a substantial component of low-energy states (including |3^)) refiecting this one-collision 
scenario. In other words, P(3^|64'^) is non- negligible because a single collision with the 
piston scatters the particle from the high-energy state 164"^) to the low-energy state |3^). 
The same argument explains, quantitatively, why the right peak occurs at ^ 320 m 



Fig-E 



Alternatively, we can use Eq. Mto rewrite P(3^|64^) as P(64^|3^), which is the prob- 
ability to end in state 164"^), starting from state |3'^) when compressing at piston speed 
V = 100. Here we imagine a particle initially moving with speed 



1^1 = ^2Ef=3 ^ 5. (26) 

As the piston moves from B = 2.0 to A = 1.0, the particle might suffer a single collision 
with the piston, imparting a leftward velocity Am ~ —2v = —200 to the particle. Thus 
for the initial state |$(0)) = |3-^) we expect the final state |$(t)) to be a superposition of 
low-energy states (corresponding to no collisions) and high-energy states near |64"^) (one 



collision). This is indeed the spectrum seen in Fig. 2(b) This interpretation suggests that 
is equal to the probability that the particle suffers a collision with the piston during the 
compression process, and is the probability it avoids a collision, when starting from state 
\n^). Semiclassically and in the limit v — oo, the probability to avoid a collision during 
compression is just the probability to find the particle in the region < x < A at time t = 
(when the piston is at location B), which leads to 

hm s^'^' = 4 = -• (27) 
v^oo Br 



The superscript "sc" emphasizes that this is a semiclassical approximation. Eq. 27 agrees 



with the term 1/r in the expression appearing in Eq. 11 (just before the inequality); the 
oscillatory term there, sin(27m/r)/27rn, is quantum-mechanical in origin. 

In either case - expansion or compression - the dynamic component is associated semi- 



classically with a collision between the particle and the piston. We conclude that Eq. 11 
underestimates s„ because it neglects the contribution due to a particle-piston collision. 
These considerations relate to the ordering of limits. Fig. [2] suggests that 

lim P{n^\m^) = 0{n^\m'^), (28) 
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for any fixed initial state Im"^). Now, Eq. 11 implicitly contains a double limit, namely, 

K 

lim s„ = lim lim P(n^|m^). (29) 

m=l 

If we take the limit K oo first (with v fixed), then both the static and dynamic components 
and are included in the sum, and the right side of Eq. 29 sums to unity (Eq. 21): 

K 

lim lim VP(n^|m^) = l. (30) 

m=l 

However, if we reverse the ordering of limits and first take f — ?■ oo (with K fixed), then 
the dynamic component gets pushed beyond the value of K, and only the static component 
contributes: 

lim lim y P(n^|m^) = lim V 0(n^|m^) = - - "^^(^^^/^) . (31) 

m=l m=l 

The physical interpretation should be clear. For any fixed piston speed v, the sudden 
approximation breaks down if m^n/A > v; for such initial states the evolving wavefunction 
catches up with the moving piston. Therefore if we sum over all initial states at fixed v, then 
this sum necessarily includes states that violate the sudden approximation. Conversely, the 



use of the sudden approximation in Eq. [TT] is equivalent to imposing a cutoff K on the sum 
over initial states: the effect of this cutoff is to exclude those states that give rise to the 
dynamic component, s^. 

The result appearing in Eq. [31] is the same as that obtained for the process of sudden 
expansion into a vacuum, in which the length of the box increases instantaneously from A to 
B. This case, considered exphcitly by TM (see Eq. 23 of Ref. [12]) and for the classical piston 
by Sung [13] , highlights the importance of the ordering of limits for the validity of Eq. [l} 
This issue is discussed in detail by Presse and Silbey [15j. See also Kurchan's lectures [20] 
for an alternative analysis of the sudden expansion process. 

While our arguments establish that Eq. [T]is valid for any piston expansion speed v, they 
also imply that for large v, transitions {n^\U\m'^) from high-lying initial energy eigenstates 



make a large contribution to s„ and ultimately to {e~^^) (Eq. [?]). When the energies of 
such high- lying states are much greater than then the probability to sample these states 
from the initial canonical distribution, P^'^ oc e"^^*", becomes exceedingly small. In this 
case, even though Eq. [T] is valid, the number of realizations required to confirm its validity is 
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FIG. 3: The probability to generate a realization from \m^) to |3^) is plotted, for the same 
parameters as in Fig. [2| and taking /3 = 0.01. 



prohibitively large. Fig. [s] illustrates this point by displaying the product P^'"^ P(n^\m'^), 
that is the net probability to generate a realization with initial and final states \m^) and 
|n^), respectively, setting /3 = 0.01 and n = 3. Comparing Figs. [2] and [sj we see that 
although realizations that correspond to the dynamic component represent an important 
contribution to s„, the probability to observe these realizations is vanishingly small. This 
conclusion is mirrored in the classical version of this expanding piston [HI [16], where a 
substantial contribution to (e~^^) arises from single-collision events, in which the particle 
loses energy as it strikes the rapidly receding piston. If Mv"^ ^ then many realizations 
of the process are needed in order to stand a decent chance of sampling initial conditions in 
which the particle is moving sufficiently fast to collide with the piston. By analogy with the 
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classical calculations of Ref . [HI [16] , we expect that the number of realizations needed for 
the convergence of the exponential average in Eq. [l] scales like exp(/3Mf^), for large v. 

We note in passing that in Figs. 2(b) and 2(c) the right peak itself exhibits a double-peak 
structure. This too has a semiclassical interpretation, which is easiest to explain in terms 
of the compression process. At t = in the state |3^), the particle is moving with speed 
|m| ~ 5 (Eq. 26). Its speed after a collision with the leftward-moving piston is greater if 
the particle was moving toward the piston just before the collision {u ~ -|-5) than if it was 
moving away from the piston {u ~ —5). A back-of-the envelope calculation suggests that this 
difference splits the right peak into two sub-components separated by Am = 2A\u\/t^ ~ 3, 



in agreement with what we see in Figs. 2(b) and 2(c' 



Finally, since this model provides a useful pedagogical illustration of a quantum nonequi- 
librium process (see also Refs. [2T1 - [2^ ). we briefly discuss the work distribution ppiW) for 



the expanding quantum piston (see Eq. 22), plotted in Fig. |4]for various piston speeds. In 
the limit f — )■ 0, the quantum adiabatic theorem gives us P{n^\m^) — )■ 5mn- Thus the work 



distribution in Fig. 4(a) reflects the initial thermal energy distribution: the largest peak 
corresponds to the situation in which the system begins and ends in the ground state, the 
next largest corresponds to the first excited state, and so on. In the opposite limit of large 
f, Pf{W) approaches an asymptotic distribution, obtained by replacing P{n^\m"^) with 



its static component 0{n \m ) in Eq. 22, (However, the dynamic component, which gets 
pushed off to infinity as discussed earlier, remains essential for the the validity of Eq. [l}) 
There are two uniquely quantal features of the distributions shown in Fig. |4j First, for v >2 
we can clearly see a nonzero probability to obtain a positive value of work. This is forbid- 
den in the classical the particle loses energy each time it collides with the piston. 
Second, for the classical expanding piston the probability to obtain W = Q approaches unity 
as f — )■ oo, whereas for the quantum piston with A = 1.0 and B = 2.0 this probability 



approaches 1/2, as illustrated by the peak at = in Fig. 4(f) Finally, although it might 
not be obvious from Fig. |4| the average work performed in the limit f — > oo is zero for the 
quantum piston [25j, just as it is for the classical piston. 

To conclude, we have used exact solutions of the time-dependent Schrodinger equation 
to study the validity of nonequilibrium work relations (Eqs. [l| [2]) for the quantum piston, 
focusing on the limit of a rapidly expanding piston, v oo. Our investigation was motivated 
by Teifel and Mahler's study [12j, which highlighted the subtleties that arise when the 
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FIG. 4: Work distribution for the expanding quantum piston. Here A = 1, B = 2, /3 = 0.01, 
and the piston speed ranges from w = 0.1 to u = 1000. The free energy difference is AF = 
-f3-^ ln{B/A) ^ -30.10. The left tail of the distributions in the region W < —100 is not shown. 
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system's Hilbert space changes due to the motion of hard boundaries. As in the classical 
case, we found that both Eqs. [T] and |2] remain valid for any finite piston speed, but the 
convergence of (e-^'^) to e-^^^ requires a sum over all possible realizations. In particular, 
when V ^ important contributions arise from those rare realizations in which the 

particle begins with a sufficiently high energy to collide with the piston. These realizations 
show up as the dynamic component (the right peak) in Fig. |2] Although we have considered 
only the one-dimensional quantum piston, we speculate that similar conclusions will apply 
to more complicated quantum systems involving moving hard boundaries, for which exact 
solutions of the Schrodinger equation are unavailable. 
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Appendix A 



Eq. 17 gives the following expression for the transition probability from \m^) to \n^) 
during the expansion process: 
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For the contraction process, the transition probability from \n^) to \m^) is obtained from 
this result by making the replacements m ^ n, A ^ B, and v — )■ —v. 
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Comparing these expressions, it is straightforward to verify that they are equal: 



P(n^|m^) = P(m^|n^) 



Ai„B\ 



(A3) 
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